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ABSTRACT 

The yield of Earth-like planets will likely be a primary science metric for 
fntnre space-based missions that will drive telescope apertnre size. Maximizing 
the exoEarth candidate yield is therefore critical to minimizing the reqnired aper¬ 
tnre. Here we describe a method for exoEarth candidate yield maximization that 
simnltaneonsly optimizes, for the hrst time, the targets chosen for observation, 
the nnmber of visits to each target, the delay time between visits, and the expo- 
snre time of every observation. This code calculates both the detection time and 
multi-wavelength spectral characterization time required for planets. We also 
rehne the astrophysical assumptions used as inputs to these calculations, rely¬ 
ing on published estimates of planetary occurrence rates as well as theoretical 
and observational constraints on terrestrial planet sizes and classical habitable 
zones. Given these astrophysical assumptions, optimistic telescope and instru¬ 
ment assumptions, and our new completeness code that produces the highest 
yields to date, we suggest lower limits on the aperture size required to detect and 
characterize a statistically-motivated sample of exoEarths. 

Subject headings: telescopes — methods: numerical — planetary systems 


1. Introduction 


The direct detection of Earth-like planets will be a primary science goal for future 
space-based flagship missions and will likely drive telescope and instrument design, includ¬ 
ing aperture size. Many studies have estimated the number of Earth-sized pl anets residing 


in the habitable zone (HZ) that could be detected with future missions (e.g., iBrowru 12005 
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Brown fc Soumm^boiO : Savranskv et alj[2010 : Turnbull et al. 2012 : Stark et all2014 ). How¬ 
ever, these e stimated yie l ds var y by a factor of a few, even when adopting similar mission 


capabilities. IStark et ahl fl2014j) showed this variation is due in part to differences in the 


exposure time distribution among candidate stars, and provided an Altruistic Yield Opti¬ 
mization (AYO) method to maximize exoEarth candidate yield by optimizing the exposure 
time distribution. Remaining yield discrepancies are the result of differing astrophysical as¬ 
sumptions, observational requirements, and the ability of the yield calculation code to fully 
simulate the execution of a mission. 


Here we build on the work of IStark et ahl (120141) to address these issues and update our 
calculations to include wavelength-dependent flux calculations, multiple visits per star, and 
spectral characterization time estimates. We hrst address the astrophysical assumptions 
made when calculating exoEarth yield and provide suggested values derived from current 
published results. In Section 12.21 we discuss our scientihc motivation, as well as the obser¬ 
vational, telescope, and instrument assumptions made. In Section [3] we provide new yield 
optimization methods to incorporate multiple visits to each star, further maximizing ex¬ 
oEarth candidate yield. We quantify the impact of these new methods and use them to 
estimate exoEarth candidate yield as a function of mission parameters in Section 01 Finally, 
in Section [5] we discuss our results, dehne a minimum scientihcally compelling yield goal, 
and use this yield goal to provide rough constraints on the required aperture size. 


2. Assumptions & Caveats 
2.1. Astrophysical Assumptions 


The exoEarth candidate yield for a direct imaging mission is dependent on several 
assumed astrophysical quantities. Most notably among these are the assumed size and optical 
properties for every Earth-sized planet residing in the habitable zone, the exozodiacal dust 
cloud surface brightness, and the number of habitable zone Earth-sized planets per star 
Here we briefly summarize our assumptions in each of these areas. 


While certainly not all exoEarth candidates will be similar to Earth in size, composition, 
or albedo, we have no observational constraints on this variety. Thus, most studies to date 
have simply adopted an Earth-twin as the dehnition of an exoEarth ca ndidate, tho u gh th e 
geometric albedo of this Earth-twin varies from study to study. As in IStark et ahl (120141 ). 
we continue to adopt an Earth-twin with a V band geometric albedo of 0.2 as our dehnition 
of an exoEarth candidate. Table [T] lists the parameters of our Earth-twin. 


We assume the same dehnition of 1 “zodi” of exozodiacal dust as IStark et ahl (120141) : 
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a fixed optical depth that produces a V band surface brightness of 22 mag arcsec”^ around 
a Solar twin (see Appendix C.2 in IStark et ahl f)20141) i. but also calculate the wavelength- 
dependent flux. W e also adopt the same local zodi flux calculations (see Appendix B in 
Stark et ahl (120141) 1. but include the wavelength dependence of the zodiacal cloud. In Section 


13.21 we discuss our calculation of the wavelength dependence of both of these sources. 

The most significant change to our astrophysica l assumptio ns relates to ou r definition of 
the clas sical ha b itable zone iHZl. lStark et ahl (120141) adopted the iBrownl (120051) HZ definition 
(BHZ). iBrownI (120051) distributed orbital semi-major axes from 0.7 ~ 1.5 AU for a Sun-like 
star, scaling each se mi-maior axis by -v/L^/L m. Detailed estimates of the classical HZ are 
slightly more distant. iKopparapu et ahl (120131) presents several limits to the inner and outer 
edges of the classical HZ, with an “optimistic HZ” (OKHZ) ranging from 0. 75 - 1.77 AU and a 
“pessi mistic HZ” (PKHZ) ranging from 0.99 - 1.67 AU for a Sun-like star. IKopparapu et ah 


( 2013 ) shows that the classical HZ deviates by a few percent from the \JL^/Lq scaling 
relation, though we continue use of this simple scaling relation out of expediency. 


Additionally, iBrownI (120051 ) distributed orbits uniformly in se mi-mai o r axis and eccen¬ 
tricity from 0-0.35. As a result, the exoEarth candidates in IBrownI (120051) can range 
in circumstellar distance from roughly 0.5 - 2.0 AU for a Sun-like star. We adopt semi¬ 
major axes distributed uniformly i n logarithmic space, as supported by radial velocity sur¬ 


veys a nd Kepler observations (e.g.. lCumming et al.ll2008l: iHoward et al.ll2012l:IPetigura et ah 


2013bl). Th o ugh w e have little direct evidence of the eccentricities of HZ Earth-sized planets, 
Kane et ahl (120121 ) showed a significant trend in the average eccentricity derived from Ke¬ 
pler data. By comparing deviat ions in the transit durations as a function of planet radius to 
model orbits, iKane et ahl (120121) showed that planets smaller than Neptune have increasingly 
circular orbits, likely due to the fact that smaller planets tend to be in multi-planet systems. 
We therefore adopt circular orbits. 


The assumed orbital distribution of planets can lead to differences in expected yield. 
For illustrative purposes. Figure [1] shows example yields as functions of aperture size for six 
different assumptions about the HZ orbits, assuming a hxed value of r;®. These different 
definitions result in a ~ 20% spread in yield. 
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Fig. 1.— Sample yields as functions of telescope diameter for six different distributions of HZ 
orbits, assuming a fixed valued of r/®. The assumed orbital distribution alone can produce a 
~ 20 % spread in yield. 


However, each HZ dehnition should have its own corresponding value for 77 ®. Recent 
estimates of 77 ® for Solar type stars suggest 77 ® gs 0.2 for planets with r a dii 1-2 Rm receivin g 
0.25 - 4 times t he Ea rth’s incident Solar flux f Petigura et al.ll2013al: ISilburt et ahl 12014 ). 


Petigura et ahl fj2013al) also found that 77 ® scales with the area of the box with dimensions 


log (Rpianet,max/-^planet,min) ^ log ( Uonter etijyR /Qinner ^r p-e)- T huS, tO Calculate 8 ^alue for 77® for 
Sun-like stars using the results of IPetigura et ahl fl2013al) and lSilburt et ahl (120141) . we must 
also define the range of planet radii that we consider “Earth-like.” 

Estimates of the bulk densities of planets as a function o f radius show a transition at 


1.5 R® (e.g., IWeiss fc Marcvl 120141 IWolfgang fc LopezI 120141) . Planets smaller than ~ 1.5 


R® are likely rocky in composition while larger planets are likely enveloped by a H/He 
atmosphere. We therefore adopt 1.5 R® as the maximum exoEarth candidate radius. 

Although there is no firm lower-limit o n the size of a pote n tially habitable world, loose 
guidelines can be derived from the work of lZahnle fc Catling! (120131) . This work proposes 
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a lower limit on the mass of habitable planets given the inpnt of high energy radiation to 
the planet and theoretical atmospheric loss processes. For a planet that receives the same 
flnx from its star that Earth receives from the Snn, the lower limit is ~ 0.6 assnming an 
Earth-like bnlk density. 


However, we also reqnire that onr Earth-twin assnmption adopted for yield calculations 
be representative of the typical exoEarth candidate radius. Thus, we demand that the median 
planet radius be equal to 1 R^. This results in our adop ted minimum radius of 0 .66 i? 0 , 
roughly consistent with the minimum radius suggested by Zahnle fc Catling (2013 ). Given 
this range of planet radii, we scaled from the results of iPetigura et al.l fj2013al) for each 
of the three classical HZ dehnitions discussed, maintaining a constant fractional uncertainty. 
The results are presented in Table [2J 


Figure [2] shows the same sample yields shown in Figure [H but scaled by the //q values 
listed in Table [2l The spread in the yield between the BHZ and the OKHZ definitions 
decreases to ~ 10%. The yield associated with the PKHZ is ~ 40% smaller. 



Fig. 2.— Sample yields as functions of telescope diameter for six different distributions of 
HZ orbits, assuming the self-consistent values of listed in Table |2J 
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Table 1. Baseline Astrophysical Parameters 


Parameter 

Value 

Description 

»?© 

0.1 

Fraction of Sun-like stars with an exoEarth candidate 

Rp 

1 

Planet radius 

a 

[0.75,1.77] AU* 

Semi-major axis (uniform in log a) 

e 

0 

Eccentricity (circular orbits) 

cos i 

[-1.1] 

Cosine of inclination (uniform distribution) 

U) 

[0, 2n) 

Argument of pericenter (uniform distribution) 

M 

[0, 27r) 

Mean anomaly (uniform distribution) 

# 

Lambertian 

Phase function 


0.2 

Geometric albedo of planet at 0.55 and 1 /rm 

z 

23 mag arcsec“^^ 

V band surface brightness of zodiacal light 

X 

22 mag arcsec“^^ 

V band surface brightness of 1 zodi of exozodiacal dust 

n 

3 

Number of zodis for all stars 


*a given for a solar twin. The habitable zone is scaled by L* /Lq after calculating 

projected separation Sp. 

^Varies with ecliptic latitude—see Appendix B in IStark et al. I JSoIi) . 

^For Solar twin. Varies with spectral type—^see Appendix C in IStark et al. li^. 


Table 2. t]® for Planet Radii 0.66 - 1.5 


HZ Definition 

Acronym 

a- 1 

limner 

(AU) 

Clouter 

(AU) 


Brown (2005) 

BHZ 

0.7 

1.5 

0.14 ±0.05 

Optimistic Kopparapu et al. (2013) 

OKHZ 

0.75 

1.77 

0.16 ±0.06 

Pessimistic KoDoaraou et al. (2013) 

PKHZ 

0.99 

1.67 

0.10 ±0.04 


^For a Sun-like star. 


^Based on results and scaling relations from l^tigura et al.l 
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The choice of which classical HZ to adopt is a judgement call. Designing a mission 
around a pessimistic HZ will lead to a larger required aperture and will more hnely tune the 
mission to a smaller region of orbital space. Choosing an optimistic HZ may result in more 
Earth-sized planets with a smaller aperture, but fewer truly Earth-like planets. We note that 
the la lower limit on r]^ for the OKHZ is equal to the expected value of 77 ® for the PKHZ. 
For the purposes of this paper, we choose the OKHZ definition with planets on circular orbits 
distributed uniformly in the logarithm of semi-major axis, and use the conservative la lower 
limit on of 0.1. 


We based our 77 ® assumption on the work of IPetigura et ahl fj2013al) . Howev er, other re¬ 


cent planet occurrence rate es 


2013 


Foreman-Mackey et al 


imat es suggest somewhat smaller values of 77 ® fe. g..lDong &: Zhu 


2014). Ev en when adopting the same catalog as iPetigura et al. 


f 2013a '). Foreman-Mackey et ahl f 20l4 derives an 77 ® value three times smaller than that 


Petigura et al.l (120134 . We therefore do not regard our adopted value of 77 ® = 0.1 as 


of 


pessimistic. 


2.2. Telescope, Instrument, &: Observational Assumptions 


We calculate exoEarth candidate yield for a coronagraphic mission only, ignoring alter¬ 
native high-contrast imaging technologies like an external starshade. We calculate yield via 
Monte Carlo simulations of completeness. We ignore scheduling constraints and the effects 
of any real-time decisions that may be made during the course of the mission. 

We also ignore the impacts of overheads associated with telescope slews and wavefront 
control, which are architecture-dependent. The overheads for a coronagraph may be as long 
as the observations themselves. Thus, in reality the required mission lifetime may be twice 
as long as the assumed total exposure time of the simulations in this work. Additionally, 
pointing overheads for the telescope may signihcantly reduce the yields of the largest aper¬ 
tures we simulate, for which the slew time may be comparable to some of the exposure times. 
We leave the addition of overheads to future work. 


We assume a mission that has the primary goal of detecting an Earth-like planet. The 
mission must therefore be capable of remotely characterizing any exoEarth candidates to 
determine whether or not they are Earth-like. There is no agreed-upon dehnition of “Earth- 
like,” but one commonly accepted minimum requisite is liquid surface water, potentially 
resulting in detectable amounts of atmospheric water vapor. The water va por feature that is 
comm only viewed as the most readily detectable occurs near 0.95 pm fe.g.. lBrandt fc Spiegel 


2014 . We therefore adopt the search for atmospheric water vapor at 1 pm as an observational 



































requirement. While we continue to assume V band detections of exoEarth candidates, only 
the exoEarth candidates directly observable at a wavelength of 1 fim count toward the total 
yield of the mission. 


We have updated the code of IStark et ahl (120141) to include the time required for spectral 
characterization of all detected exoEarth candidates at 1 pm. The spectral characterization 
time counts against the total exposure time budget of the mission, i.e. a higher detection rate 
results in more time spent characterizing. For the baseline mission, we require an i? = 50 
spectrum with S/N = 5 per spectral resolution element at a wavelength of 1 pm for every 
detect a.ble exoEarth candidate , which provides a 50% probability of a 3a detection of water 
vapor flBrandt fc Spiegel! 120141) . We implicitly assume that all non-exoEarth candidates (all 
planets outside of the HZ and all planets within the HZ that are not Earth-sized) can be 
distinguished from exoEarth candidates without taking i? = 50 spectra, either by brightness, 
color information, or lower resolution spectra that require significantly less exposure time. 


The search for water defines our minimum characterization requirement. Searching 
for additional biosignatures may require additional exposure time if multi-band observations 
cannot be made in parallel using multiple coronagraphs. However, this additional time could 
be reduced by applying a “ladder of biosignatures,” in which only the exoEarth candidates 
that exhibit detectable traces of water are followed up for additional characterization. 


Unlike IStark et ahl (120141) . we include the effects of detector noise. We assume noiseless 
detectors for the baseline mission and address the impact of detector noise in Section We 
also approximately include coronagraphic broadening of the planet’s PSF via the broadening 
factor b, which linearly scales the area of the planet’s Airy pattern PSF. We assume a 
diffraction-limited planet PSF {b = 1) for the baseline mission. 


Finally, we allow stars to be observed multiple times, i.e. revisited, to increase the 
chances of detecting exoEarth candidates. We place no restrictions on the exposure time 
of each observation. We do not include any revisits that may be required to establish the 
orbit of a potentially habitable planet. Depending on the observation plan implemented, 
orbit-constraining revisits may be rarely required, and we expect their impact on the overall 
mission yield to be small given that most stars are already assumed to be observed multiple 
times for the purposes of increasing the yield. 

Table [3] lists the baseline telescope and instrument parameters. We assume a telescope 
of diameter D paired with two coronagraphs, each with its own detector, whose parameters 
are chosen optimistically to represent future performance. The first coronagraph is designed 
for broadband detection at Ad = 0.55 pm. We will refer to this hereafter as the “detection 
coronagraph” and subscript its dehning parameters with a ‘d.’ Because p® = 0.1 implies a 
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< 10% chance of detecting an exoEarth candidate during any observation, many observations 
must be made, motivating the need for an optimistic raw contrast of 10“^° to reduce detection 
exposure times. 


The second coronagraph is designed for spectral characterization at Ac = 1 /xm. We 
will refer to this hereafter as the “characterization coronagraph” and subscript its defining 
parameters with a ‘c.’ Because the inner working angle of a coronagraph roughly scales with 
\/D and we require the detection of water for every exoEarth candidate, the inner working 
angle of our 1 /rm characterization coronagraph defines which planets are observable. We 
thus choose the same physical IWA for both coronagraphs, with an optimistic IWAc = 2Xc/D 
and a relaxed IWAd = 3.6Xd/D. Because the two inner working angles are the same physical 
size, we will refer to them simply with the singular parameter IWA. 


Realistically, the optimistic IWAc = 2Xc/D may have to come at the expense of con¬ 
trast. Because we assume spectral characterization only occurs for true exoEarth candidates, 
spectral characterization will occur on relatively few targets. Thus, the characterization time 
may not greatly impact the overall yield. In light of this we relax the baseline raw contrast 
requirement at 1 /xm to Cc = 5x10“^*^. Both coronagraphs are also assumed to have the same 
physical OWA as well as the same absolute noise floor. In summary, we adopt a coronagraph 
for 1 fim spectral characterization that is optimized for small IWA, and a coronagraph for V 
band detection that is optimized for high contrast/short exposure times. 


We assume a total exposure time of 1 year out of a 5 year mission (though we do not 
explicitly budget for overheads, 1 additional year may be required for telescope and corona¬ 
graph overheads). We assume the faintest point source that can be detected at the threshold 
S/N is 2.5 magnitudes fainter than the raw detection contrast Cd, i-e. post-processing can 
reveal point sources with relative flux equal to one tenth that of the raw detection contrast, 
and one fiftieth that of the raw spectral characterization contrast for the baseline values of 
Cd and Cc- Second, we relax the threshold S/N for detection from a value of 10 to 7, a choice 
we explain in Section 14.31 


We also make several minor changes to the root target list. In IStark et al.l (120141) . we 
selected targets within 50 pc from the Hipparcos catalog, removing only those targets with 
known companions within 10". Here, we also remove stars categorized as luminosity classes 
I, II, and III as well as stars that have no classification, to select only main sequence and 
sub-giant stars. These changes are negligible, impacting the calculated yield at the 1% level. 


We note that we did not vet the input target list in great detail. We vetted binaries 
by simply removing those stars with companion separations < 10", where companion sep¬ 
arations were taken as the minimum of the separations reported in the Hipparcos Double 
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& Multiples Catalog and the Washington Double Star Catalog. Additionally, the Hippar- 
cos catalog is missing several nearby, bright stars that would make excellent targets for 
exoEarth searches. A more thorough vetting of candidates would reference additional cat¬ 
alogs for missing stars, remove any additional bad targets from the Hipparcos input list, 
and estimate future binary configurations. However, under the assumption of = 0.1, 10 
additional/fewer high priority targets with fully observable habitable zones would change 
the yield by less than 1 exoEarth candidate, so we expect the impact of these additional 
screenings on the mission yields reported in this paper to be minimal. 


3. Methods 


The fundamentals of our exoEarth c andidate yield calculations are based on iBrown 


(120051) and described in IStark et al.l (120141 . Here we discuss updates to the code. 


3.1. Detector Noise 


Stark et ahl (12014) calculated the required exposure time assuming negligible detector 


noise. Here we explicitly include the detector noise count rate as an additional parameter to 
constrain acceptable levels. We calculate the background count rate as 


CRb 


CRb,* + CRb,zodi + CR 


^b,exozodi 


CR, 


b,detector ? 


( 1 ) 


where CRb,*, CRb,zodi, and CRb,exozodi are the count r ates of the unsupp ressed stellar light, 
local zodi, and exozodi, given by Eq. 11, 13, and 15 in IStark et ahl (12014) . respectively. The 
additional term, CRb,detector, represents the noise count rate of the detector, a parameter we 
control. We can relate read noise RN and dark count rate ^ to CRb,detector via 


CRb ,detector ~ '^pix ('^ T RN /Texpose) , 


( 2 ) 


where Upix is the number of pixels required for imaging or spectral characterization and 
^expose is the length of an individual exposure. In practice, we assume the detection and 
characterization coronagraphs have independent detectors and allow for corresponding values 
of CRb,detector,d and CRb,detector,c, respectively. The above expressions underestimate the read 
noise contribution for observation times r < Texpose, but this only occurs for observations of 
bright targets, for which the unsuppressed starlight will dominate the noise sources. 
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3.2. Spectral Characterization Time 


We calculate the required spectral characterization time, tc, for every synthetic exoEarth 
around every star using the same basic exposure time expressions given in Section 3 of 
Stark et al.l (120141) . However, we use the updated expression for CRb given above, substitute 
S/N(. for S/N, and substitute AAc = Xc/Rc for AA. 

We also adjust the flux of every source to the desired spectral characterization wave¬ 
length of 1 fim. We assume a wavelength-independent geometric albedo for Earth, such that 
the planet’s flux is linearly proportional to the stellar spectrum. This assumption is valid 
for this work , as the Earth’s conti nuum geometric albedo is approximately 0.2 between 0.5 
and 1.2 /im flRobinson et al.l 1201 ll) . We calculate the flux of each star at 1 /im by linearly 
extrapolating the stellar fluxes from the V and I magnitudes given in the Hipparcos catalog. 
Although this extrapolation to 1 /rm is a rough estimate, it suffices for our calculations, as 
the yield is a weak function of total exposure time and thus neg ligibly impacted by modest 
changes to the spectral characterization time fjStark et al.ll2014j) . 


For the l o cal zo di, we calculate the wavelength dependence by interpolating Table 19 of 
Leinert et al.l (1199811 . We calculate the 1 fim brightness of the exozodi and flux of the planet 
around each star by multiplying the V band flux of each source by F*(l/xm)/F*(0.55/im), 
where F*(l /rm) and F*(0.55 pm) are the stellar fluxes at 1 pm and 0.55 pm, respectively. 


Given that the spectral characterization time will count against the total exposure time 
budget, such a mission cannot afford to take Rc = 50, S/Nc = 5 spectra of every system. 
However, it would be risky to wait to take spectra, as the planet could disappear interior 
to the IWA, move into crescent phase, or become confused with another object. Thus, the 
choice to spectrally characterize a planet must be made in near-real time. 


The most robust treatment of spectral characterization therefore requires some sort 
of Monte Carlo mission execution simulator, run many times over to statistically bound 
the impact of spectral characterization time. Such a simulation would be computationally 
expensive compared to our yield estimation code and is more appropriate for modeling a 
small number of mission designs. In lieu of robust mission execution simulations, we include 
the impact of spectral characterization time in a probabilistic sense that approximates a 
large suite of Monte Carlo simulations. 


To approximate the required spectral characterization time of the mission, we require 
a spectrum for every detected exoEarth candidate. We start by calculating tc for every 
synthetic exoEarth simulated around each star. We calculate the probabilistic spectral char- 
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acterization time for the star as 


V(B (^c)) 


( 3 ) 


where Ci is the completeness of the star, the product is the probability of detecting 
an exoEarth candidate around the star, and (fc) is the average spectral characterization 
time of all detectable synthetic exoEarths around the star. Because 


_ ^ 2 ,detectable 

i : 

^simulated 


( 4 ) 


where n,,detectable is the number of detectable synthetic exoEarths around the star and 
^simulated IS the total number of simulated synthetic exoEarths around each star, and 


detectable 


(tc) — ^ 


i=i 


^ 2 ,detectable 


we can rewrite Equation [3] as 




V® 


detectable 


^simulated 


i=i 


‘'3,Cl 


( 5 ) 


( 6 ) 


where the quantity r/®/^simulated is the probability of observing an individual simulated syn¬ 
thetic exoEarth and the summation is simply the total of the spectral characterization times 
for all detectable synthetic exoEarths. We note that the spectral characterization time counts 
against the mission’s total exposure time budget; missions that detect more planets must 
devote the required spectral characterization time to them. 


3.3. Optimized Revisits 


The majority of completeness calc ulations thus far hav e assu med a single visit to each 
star, the most notable exception being iBrown fc Soiimmerl (120101 ). who include revisits via 
an analytic approximation. The single-visit assumption is commonly made for two reasons. 
First, the initial visit to a target star is the most productive visit, since no portion of 
orbital space has been ruled out. Strictly speaking, any revisit ca n only achieve a maximum 
completeness of one minus the previous visits’ total completeness (IBrown fc Sonmmerll2010l ). 
However, an additional visit to a high priority star may prove more fruitful than the hrst 
visit to a lower priority star. 


Second, revisits require computationally intensive calculations. The orbit of every un¬ 
observed synthetic planet must be advanced by the delay time between revisits. At, prior 
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to calculating the revisit completeness. Further, one would ideally optimiz e the delay time 
At to make the revisit as efficient as possible. iBrown fc Soummerl fj2010l) estimated that 
the revisit completeness of a single star calculated for a single delay time requires on the 
order of a few seconds on a single processor; optimizing the delay time f or hundreds of stars, 
each w ith multiple visits, could take hundreds of hours. In light of this, IBrown fc Soummer 
(I 2 OIOI) opted for a non-optimized approximation to the revisit completeness. 


We have overcome the numerical run time limitations of revisits, enabling optimization. 
We accomplished this by indexing the orbits of every simulated planet in mean anomaly, 
avoiding on-the-fly orbital calculations, and prioritizing revisit calculations to avoid unnec¬ 
essary computations. We can simultaneously optimize the number of revisits to each star, 
the delay time between revisits, and the exposure time of every visit to maximize mission 
yield in a reasonable run time. For example, simulating 10® synthetic exoEarths around 
every star, our code can optimize up to 10 visits to thousands of potential target stars in 
~ 5 minutes on a single processor. 


To calculate revisits, we divide each synthetic exoEarth orbit into hm steps in mean 
anomaly. Each of these steps in mean anomaly counts as a potentially observable planet. 
For example, to simulate 10® synthetic exoEarths, we sample 10^ orbits 100 times each, i.e. 

^planets ~ ^orbits X Hm ■ 

For each star, we calculate the projected separations and required exposure time of all 
Upianets plauets ffi advauce. At any future time, we can then calculate a given planet’s evolved 
projected separation and exposure time by referencing the mean anomaly index that most 
closely corresponds to a delay time of At via 


M{At) = n At — M(t = 0), 


(7) 


where n is the mean motion. To calculate the mean motion of every orbit, we make the 
approximation (M*/Mq) where M* and L* are the stellar mass and luminos¬ 

ity, respectively. We note that the validity of this scaling does not impact any of the results 
presented here, and will only become important when determining a real world observing 
schedule. 


m 


We then calculat e the hrst visit completeness as a function of exposure time as described 


Stark et al.l (120141 . Once the completeness, Ci, and exposure time, ri, for the first visit 


is chosen by optimizing the first visit exposure time, we follow these steps to optimize the 


visit to a given star: 


1. Mark all planets with IWA < Sp(Afj_i) < OWA and Tp{Ati-i) < rj_i as having been 
observed, where Sp and Tp are the projected separation and exposure time for planet p. 










14 


2 . Advance all unobserved planets by Atj. 

3. Select all advanced, unobserved planets with IWA < Sp(Atj) < OWA and sort them 
by Tp{Ati) 

4. Calculate Ci(rj) and dCi/dTi. 

5. Optimize r* (see below). 

6 . Repeat steps 2-5 for multiple values of Af*. Choose Afj such that Cj/rj is maximized, 
thereby optimizing both Afj and r*. 


A range of visit delay times Atj is chosen that allows planets to evolve in their orbits 
without losing phase coherence. Ideal values of Afj should correspond to the time for un¬ 
observed crescent phase planets to evolve to gibbous phase or planets behind the IWA to 
become visible. Thus, the optimal value of Afj is a fraction of the median orbital period. In 
practice, we hnd 20 uniformly spaced values of At* from 0 to the median orbital period well 
resolves the optimum delay time. 


The procedure above requires optim izing the exposure time of every visit, including the 


hrst. The AYO algorithm introduced in IStark et ah 


(12014J ) provides a numerically efficient 


method for optimizing the exposure time of every star in the single-visit case by redistributing 
“packets” of exposure time among stars. However, because the completeness of each revisit is 
a function of the planets left unobserved from previous visits, redistributing a single “packet” 
of exposure time from one star to another requires recalculating Cj(rj) for all subsequent 
visits. This in principle could be done, but would likely be computationally inefficient, since 
calculating Ci{Ti) dominates the run time of the code. 


Instead we opt for the alternative “equal-slope” method introduced by iHunvadi et ah 


(120071) . Briefly, the equal-slope method works under the assumption that redistributing 
“packets” of exposure time to stars with the largest dC/dr will result in all selected stars 
having equal values of dC/dr. One must therefore guess at the optimal value of dC'/dr, 
solve for the r, at which dCi/dri = dC/dr^-aess for each star, then down-select to the set of 
observations that maximize the total completeness and £t within the total exposure time of 
the mission. This is repeated until the optimized value of dC/dr is found, a process that 
typically converges in ~ 10 iterations. 
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Fig. 3.— Top-left: Five sample orbits around HIP 27072 seen edge-on by an 8 m telescope 
at a: = -foo, resolved into 100 mean anomaly steps. IWA and OWA shown in gray. Planets 
detectable during the hrst visit, primarily in gibbous phase, are shown in red. Top-right: 
First visit planets evolved, detectable second visit planets shown in orange. Bottom-left: 
Planets continue to evolve as third visit is made (blue). Bottom-right: Detectable planets 
after four visits. 


Figure [3] illustrates the revisit method at work for our baseline mission observing HIP 
27072, an F6V star at 9 pc with a ~1000 AU K2V companion. For this illustration, we have 
used the BHZ with eccentric orbits allowed. The telescope is assumed to be at a; = -|-oo, 
pointed to the left. The gray shaded regions mark the IWA and OWA, interior and exterior 
to which planets cannot be detected, respectively. We select hve sample orbits, all seen 
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edge-on by the telescope, that are resolved to 100 steps in mean anomaly. Planets with 
positive values of x appear to be crescent phase, while those with negative values appear to 
be in gibbous phase. 

In the top-left panel, planets that could be observed during the hrst visit, mostly in 
gibbous phase, are marked with red dots. Between the IWA and OWA, the transition 
between red and open circles traces a contour of constant planet brightness. In the top-right 
panel, the planets have evolved in their orbits for a time equal to At 2 , at which time the 
most unobserved planets appear in gibbous phase, triggering the second visit (orange). The 
full simulatioil^ ends in the bottom right panel, which shows a total of four visits calculated 
for this star. 

Figure 0] shows the calculated bene£t-to-cost ratio for this set of observations. The 
hrst visit completeness Ci = 0.46. The second through fourth visit beneht-to cost ratio are 
plotted as functions of delay time. The code selects the delay time at which the beneht- 
to-cost ratio is maximized, to achieve the most completeness for as little exposure time as 
possible. Each peak is labeled with the value of completeness obtained at that time. We 
obtain 94% completeness in four visits to this star. 


^Available at www.starkspace.coin/revisits.inp4 
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Fig. 4.— Benefit-to-cost ratio as a function of delay time for the four visits selected for HIP 
27072. Completeness attained for each visit is indicated near the peak cost-to-benefit ratio. 
Optimizing the delay and exposure time of revisits can double the benefit-to-cost ratio. 


4. Results 


4.1. The impact of optimized revisits 


Stark et ahl (120141 ) did not include revisits. Here we show the change in yield when 


we “turn on” revisits, but do not include spectral characterization time. The left panel of 
Figure [5] shows the yield as a function of aperture size with and without revisits. For this 
plot, we have ignored spectral characterization time. The dotted line shows the single visit 
yield while the solid line shows the yield when allowing up to 10 visits per star. 


Optimized revisits increase yield by 35-75%, with the largest changes occuring for 
smaller apertures. More detailed calculations show that revisits increase yield by ~ 40% 
across a broad range of IWAs, contrasts, and exozodi levels. Roughly one quarter of this 
increase specihcally comes from optimizing the delay time between revisits. As we will show 
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in Section H73l yield scales roughly as D'^. Thus, the addition of optimized revisits reduces 
the required aperture by ~ 15%. 

The right panel of Figure [5] shows the number of unique stars observed. When including 
revisits, the number of stars visited decreases, as expected since stars are visited multiple 
times. The number of observations increases, though, since nearby stars that have short 
integration times are observed multiple times. 




Fig. 5.— Comparison of ExoEarth candidate yield (left) and number of unique stars observed 
(right) as functions of aperture size for the single visit and multi-visit cases. No spectral 
characterization time is included in these calculations. 


4.2. The impact of spectral characterization time 

Figure [6] shows how the multi-visit exoEarth candidate yield changes when we “turn 
on” both revisits and our calculation of spectral characterization time. For this plot, we 
use the baseline mission parameters, for which we assume noiseless detectors. The spectral 
characterization time reduces yield by 3-7%. 
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Fig. 6.— Comparison of ExoEarth candidate yield (left) and number of unique stars observed 
(right) as functions of aperture size for the cases of no spectral characterization time and an 
Rc = 50, S/Nc = 5 spectrum obtained on every candidate exoEarth. All calculations allow 
for up to 20 visits per star. 


4.3. Sensitivity of exoEarth candidate yield to mission parameters 

Here we vary one parameter at a time to show how the yield responds to deviations 
from the baseline mission. For these results, we include optimized revisits and spectral 
characterization time using the parameters listed in Table [3l 

Figure [7] shows how the exoEarth candidate yield scales with aperture as well as total 
exposure time for the mission. The red dot represents the baseline mission parameters. Each 
plot is inset with a value for 0, which is the approximate sensitivity of the yield with respect 
to the parameter being varied, i.e. the exponent of the power law relationship between yield 
and the parameter plotted, at the location of the red dot. 
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Fig. 7.— ExoEarth candidate yield as a function of mission parameters. Yield is a strong 
function of aperture and a modest function of total exposure time/mission lifetime. 


As shown on the left, near the baseline mission parameters, yield scales roughly as D^] 
the mission yield can be doubled by increasing the aperture size by ~ 40%. This is the 
strongest dependence of any parameter that we investigated, even if one considers yield to 
be a function of D'^, which would reduce 0 by a factor of two. This strong dependence is due 
to the fact that an increase in D produces three benehts: an increase in photon collection 
rate, a decrease in PSF size (reducing the amount of background noise blended with the 
planet’s signal), and a decrease in the on-sky IWA obtained with the coronagraph. Each of 
these three benehts contribute roughly equally to the dependence of yield on D. 


As shown in the right panel of Figure [71 the yield is a modest function of total exposure 
time (i.e., mission lifetime), with yield scaling as (Sr)°'^^. A doubling of mission lifetime 
only increases yield by ~ 25%. This modest dependence is a result of assuming targets 
can be observed in priority rank order. As the mission progresses, one must observe lower 
priority targets. Equivalently, if the mission must sacrihce tot al exposure t ime, t he lowest 
priority targets are removed. This selection effect, discussed in IStark et ah! (120141) . leads to 
moderately weak dependencies for many of the mission parameters that follow. 


Figure |8] shows the yield as a function of mission parameters that are most relevant 
to the internal coronagraph instrument. The top-left panels show the yield as a function 
of the coronagraphs’ IWA and OWA. We remind the reader that we choose an IWAd = 
3.6Ad/71 for the detection coronagraph at Ad = 0.55 pm, and an IWAc = 2Xc/D for the 
characterization coronagraph at Ac = 1 pm, such that both coronagraphs have the same 
physical IWA. Compared to other instrument parameters, the yield depends quite strongly 
on the coronagraphs’ IWA near the baseline mission parameters, with yield oc IWA“°'®®. An 
increase in the coronagraph IWA from the assumed 3.6\d/D (flXc/D) to bAX^/D {3Xc/D) 
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would reduce yield by ~ 35%. 
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Fig. 8.— ExoEarth candidate yield as a function of instrument parameters. 


The top-right panels in Figure |8] show the yield as a function of the coronagraphs’ 
contrasts, which are independent quantities. Yield is a weak function of the detection coro- 
nagraph’s contrast, with yield oc An order of magnitude degradation in contrast from 

the assumed 10“^° to 10“® reduces y ield by ~ 40%. Th is is slightly weaker than the depen¬ 
dence for the baseline parameters in IStark et al.l fl2014j) , primarily because here we assume 
a noise floor 2.5 magnitudes fainter than the raw contrast, such that Earth twins are ob¬ 
servable even at contrasts of Cd = 10“®. In other words, the dependence of yield on contrast 
is weak in part because we assume the unsuppressed starlight contributes noise that can 
be overcome with longer integration times. Contrast also only impacts nearby targets; the 
background count rates are dominated by local zodi and exozodi for targets farther than ~ 7 
pc, where most planets are found. 


The yield is even less sensitive to the characterization coronagraph’s contrast. This is 
because we only require spectral characterization for those planets that are truly exoEarth 
candidates, i.e., exoEarth candidates can be distinguished from other types of planets by 
brightness, orientation, and color information. If this assumption is invalid and spectral 
characterization must be performed for all observations, we expect the characterization coro¬ 
nagraph’s contrast to be at least as important as the detection coronagraph’s contrast. 
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The bottom-left panel in Fignre [8]shows the yield as a fnnction of the systematic noise 
floor. The systematic noise floor is plotted as the “limiting magnitnde,” or the difference in 
magnitndes between the sta r and the d i mmes t planet that can be detected at the desired 
threshold S/N. As shown in IStark et ahl (120141) . if one’s goal is to hnd Earth twins, there is 
no beneht to observing planets fainter than 26.5 magnitndes dimmer than their host star. 


The bottom-middle panels in Fignre [8] show the yield as a fnnction of thronghpnt and 
bandwidth. These two qnantities are mathematically eqnivalent in the detection exposnre 
time eqnation, so we may expect them to have identical sensitivities. However, nnder onr 
cnrrent assnmptions, the thronghpnt factors into both detection time and characterization 
time, whereas the bandwidth only factors into detection time. The dependence of yield on 
thronghpnt is ronghly eqnal to that of total exposnre time, as expected since they both con¬ 
trol the nnmber of photons collected within the mission lifetime. The priority rank selection 
effect previonsly discnssed prodnces the modest dependence on both of these parameters. 
We note that in reality, bandwidth will affect the wavefront control overheads (which we do 
not model in this paper) and thus indirectly affect the spectral characterization time. 


Finally, the lower-right panel shows the yield as a function of the planet’s PSF areal 
broadening factor, b. Again, we hnd the s ame moderately w eak dependence: yield is oc 
(l/^)0.33 py Equations 7-15 in IStark et al.l (120141) . in the background-limited 

regime in which the planet’s count rate is much less than the background count rate, exposure 
time is inversely proportional to both throughput and (1/H). A doubling of the planet’s PSF 
size will roughly reduce the yield by ~ 25%. 


Figure [9] shows the yield as a function of the detector noise count rates. We have not 
included 0 values in these plots as we assumed noiseless detectors for the baseline mission. 
For the detection coronagraph, detector noise count rates > 0.06 counts sec“^ decrease the 
yield by > 5%. For the characterization coronagraph, detector noise count rates > 0.5 counts 
sec“^ decrease the yield by > 5%. Clearly the yield depends more strongly on the noise count 
rate of the detection coronagraph. However, this does not necessarily mean that the required 
read noise and dark current of the detection coronagraph’s detector are lower. 
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Fig. 9.— ExoEarth candidate yield as a function of detector noise count rates. 


To translate the above count rates into combinations of read noise and dark current 
(RN, ^), we use Equation [2] and assume values for Texpose and Upix. We assume the length of 
an individual exposure Texpose = 1000 s to allow for the removal of cosmic ray artifacts. The 
number of pixels over which the planet’s photometric aperture is spread, Upix, depends on 
the optical layout of the instrument. 

For detection, we imagine three types of instruments that could be used: a simple 
imager (assuming spectral information is not needed for wavefront control), an imager with 
a frequency-resolving detector, and an integral held spectrograph (IFS). If the detection 
coronagraph is a simple imager or is frequency-resolving, the planet’s photometric aperture 
would be spread over ~ 4 pixels assuming imperfect centering and Nyquist sampling of the 
PSF. If the detection coronagraph is an IFS, then the planet’s photometric aperture would 
be spread over 4 lenslets, each of which is dispersed onto 6 pixels (2 in the frequency domain 
and 3 in the spatial domain) per spectral channel. Assuming R = 50, such that there are 
10 spectral channels over the 20% detection bandpass, the IFS would spread the planet’s 
photometric aperture over 240 pixels. 

For spectral characterization, we imagine two types of instruments: an imager with 
energy-resolving detector and an IFS. Similar to a broadband detection, the planet’s photo¬ 
metric aperture would be spread over 4 pixels for each spectral channel in an energy-resolving 
detector. As discussed previously, the IFS would spread the planet’s photometric aperture 
over 24 pixels per spectral channel (4 lenslets, 6 pixels per lenslet). 

Figure [10] shows the fraction of the yield realized for different instruments. The left 
panels show the impact of detector noise on yield for the V band detection coronagraph only. 
The right panels show the impact of detector noise on yield for the 1 pm characterization 
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coronagraph only. These curves were created by interpolating the yield curves from Figure [9] 
and adopting the values of Upix described above. They were calculated without any exposure 
time limitations due to stability requirements or motion of the observed planet. 


Detection 


Characterization 




Fig. 10.— ExoEarth candidate yield as a function of read noise and dark current, normalized 
to the noiseless detector case, for different detection coronagraph instruments (left) and 
characterization coronagraph instruments (right). Combinations of read noise and dark 
current to the lower-left of the solid line are nominally acceptable for the baseline mission 
parameters. All curves are sensitive to bandwidth/spectral resolving power and the required 
signal to noise ratio (see Figure [TT]). 


While these (RN, curves are roughly independent of aperture size, they are sensitive 
to bandwidth/spectral resolving power and S/N. Figure [TT] shows the 95% yield curves for 
the characterization coronagraph’s detector for several combinations of Rc and S/Nc. We 
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discuss the implications of these curves in Section [5l 

Characterization 



Fig. 11.— Curves for 95% of possible yield (see Figure [10]) for the characterization coro- 
nagraph’s detector for different combinations of Rc and S/N^. Estimates of the allowable 
detector noise parameters are sensitive to bandwidth/spectral resolving power and S/N. 


Figure WI\ shows the yield as a function of the threshold detection S/N, which shows a 
strong dependence. The threshold S/Nd is not a telescope or instrument parameter—it is a 
numerical parameter that dehnes the minimum signal to noise ratio allowable for detection. 
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In our calculations, we assume that the signal to noise ratio of an exoEarth candidate follows 
Poisson statistics; we estimate the noise by simply taking the square root of the number 
of counts. If this were truly the case, we could justifiably choose a threshold S/N^ as 
low as 5. In reality, unknown systematic uncertainties and complexities with background 
modeling/subtraction will li kely reduce the true S/N, motivat i ng previous work s to choose 


a threshold S/Nd = 10 (e.g., ISrownl 120051: iTurnbull et ahll2012l: IStark et ahl 1201411 . 



Fig. 12.— ExoEarth candidate yield as a function of threshold S/N required for detection. 
Yield strongly depends on S/Nd. 


Figure [13] shows a histogram of realized S/N for the baseline mission when adopting a 
threshold S/Nd = 10. Because our Monte Carlo calculations produce synthetic exoEarths 
covering a broad range of phase angles and circumstellar distances, the “cloud” of exoEarths 
around a given star has a broad distribution of planet brightness. Our code tunes the 
exposure time to the dimmest planet desired for detection using the threshold S/Nd. As a 
result, all other potentially detectable planets around a star, which are by definition brighter, 
will have realized S/N > S/Nd for the chosen exposure time. 
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Fig. 13.— Realized S/N for all synthetic exoEarths detectable with the baseline mission 
assnming a threshold S/Nd = 10. A threshold S/Nd = 10 is too conservative, leading to a 
median S/N = 18. 


As shown by Figure IT^ a threshold S/Nd = 10 prodnces a median S/N = 18. This 
“typical” S/N is far more conservative than necessary. In light of this, we have rednced 
the baseline threshold S/Nd for this stndy. We adopt a baseline threshold S/Nd = 7, which 
prodnces a m edian S/N = 12.5, motiv ated by the fact that S/Nd = 7 provides a low false¬ 
negative rate (iKasdin fc BraemsI 120061 ). 


4.4. Sensitivity of exoEarth candidate yield to astrophysical parameters 

Fignre flTI shows how the exoEarth yield responds to changes in astrophysical parameters. 
The left panel shows that yield is a strong fnnction of with 0 = 0.96. The yield is not 
qnite linearly proportional to 7 ^ 0 , becanse larger valnes of demand more characterization 
time at the expense of detection time. Other combinations of mission parameters reveal that 
the dependence on 7/0 can be even weaker, with cj) < 0 . 8 . 
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<Geometric albedo> Median exozodi level (zodis) 


Fig. 14.— ExoEarth candidate yield as functions of astrophysical parameters. 


The middle panel of Figure [T3] shows that yield also strongly depends on the typical 
geometric albedo of the exoEarth candidates, a parameter we likely cannot constrain prior to 
a direct imaging mission. Our baseline mission assumes a geometric albedo for all exoEarth 
candidates equal to that of the Earth. If the typical exoEarth candidate albedo is less than 
that of Earth (e.g., because many of the planets do not have cloud cover), we will miss some 
of the low-albedo exoEarth candidates, but the probability of detecting the Earth-twins 
remains the same. 


Finally, the panel on the right shows the sensitivity of the yield to changes in the 
median exozodi level. Near the baseline mission parameters (j) = —0.17, such tha t a factor of 
10 inc rease in the median exozodi level reduces yield by ~ 35%. As explained bv IStark et ah 


(120141) ■ this weak dependence is due to the priority selection effect previously discussed, 
as well as changes to the optimum observation strategy: for larger exozodi levels, planet 
searches should limit themselves to more gibbous phases. 


We note that although the right-most panel in Figure [TT] is labeled “median exozodi 
level,” we in fact set every target to the same exozodi level. Of course the universe will deal 
a distribution of exozodi levels, so it is important for us to consider the impact of such a 
distribution. To investigate this, we assigned each star in our target list a random exozodi 
level drawn from a distribution. We considered two distributions: a uniform distribution 
from 0 to 2wmp Hian , a nd a log-normal distrib ution with a median exozodi level ^median and 
a parameter of 1.25 flMennesson et al.l 120141) . We ran 15 independent instances of the yield 
maximization code, randomizing the exozodi levels each time. 


Figure US] shows the yield as a function of exozodi level for a uniform distribution (left) 
and a log-normal distribution (right). For this hgure, we adopted a 4 m telescope aperture 
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and ignored spectral characterization time to rednce the rnn time of the code. The black 
dotted line in each plot represents the yield if all stars have the same exozodi level. The 
green cnrves (and green shaded regions) represent the yield (and la deviation) if every star’s 
individnal exozodi level is known ahead of time. 




Fig. 15.— ExoEarth candidate yield as a fnnction of median exozodi level for a nniform 
distribntion of exozodi levels (left) and a log-normal distribntion (right). Updating the 
observation plan as individnal exozodi levels are measnred can mitigate any negative impacts 
of the nnknown exozodi distribntion. 

The red curve shows the yield if we know the median in advance, but no individual 
levels, and refuse to adjust the observation plan as the mission progresses. Assuming we 
know the median exozodi level in advance is a valid approximation because we should be 
able to estimate the median after a dozen or so observations. However, refusing to adjust the 
observation plan as we learn individual exozodi levels makes this a conservative estimate. 

The blue curve shows the yield if we know the median exozodi level in advance (or 
very early in the mission), and we update the observation plan as we observe each star and 
measure its exozodi level. To calculate this curve, we first ran the yield maximization code 
under the assumption that all stars have an exozodi level equal to the median exozodi level 
to obtain an initial observation plan. Then we ran the code again, allowing one star’s true 
exozodi level to be known. We locked that star’s observations into the plan, prohibiting 
them from being adjusted later, then ran the code again with one more star’s true exozodi 
level known. As each star was observed and its exozodi level was measured, we allowed the 
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code to re-optimize all unobserved stars. We assumed stars were observed in priority order 
and continued this process until the baseline total exposure time was reached. 

As shown in Figure [T^ a distribution of exozodi level does not dramatically impact the 
mission yield. Our approximation of all stars having an exozodi level equal to the median 
level is valid. As long as we allow the observation plan to adjust to individual star’s exozodi 
levels, we can avoid the negative impacts of the unknown exozodi distribution and, instead, 
obtain small increases in the exoEarth candidate yield. 


5. Discussion 


The sc a ling r elationships derived in this work are very similar to those derived in 
Stark et al.l (1201411 . Surprisingly, while revisits and spectral characterization time affect 


the absolute yield, they do not signihcantly impact the scaling relationships. 


The sensitivity parameters {(f) values) calculated in the previous section inform us of the 
approximate power law relationships between yield and each parameter when deviating from 
the baseline mission. However, they do not necessarily tell us which mission parameters will 
ultimately produce the greatest increase in yield. For example, the coronagraph IWA sensi¬ 
tivity (0 = —0.98) is signihcantly larger than the detection coronagraph’s contrast sensitivity 
(0 = —0.10). An improvement of the coronagraph IWA by just 0.5Ac/D would increase yield 
by rsj 30%. To produce an equivalent yield increase via the detection coronagraph’s contrast 
would require roughly a factor of 10 improvement. But give n current coronagrap h contrasts 
and IWAs, both of these may be equally likely/challenging fjShaklan et ahl 120131 1. 


Similarly, although the yield is weakly dependent on exozodi level {(f> = —0.17), we have 
relatively poor constraints on the median exozodi level. The current upper-limit on median 
exozodi level is 60 zodis flMennesson et al.l 120141) . In the absence of future constraints, 
we should plan a mission capable of dealing with this level of exozodiacal dust to avoid 
unnecessary risk. If we can constrain the median exozodi level to ~ 3 zodis, which may be 
possible with the Large Binocular Telescope Interferometer (LBTI), the yield would increase 
by more than a factor of two, reducing the required aperture by ~ 30%. 


The yield strongly depends on the threshold S/N^, i.e., the dehnition of a “detection.” 
In Section 01 we justihed relaxing this parameter from the commonly adopted value of 10 
to our baseline assumption of 7. This alone increased yield by 30%. Additional reductions 
may be justihable as well; reducing the threshold S/Nd to 5 would provide another 30% 
increase in yield. However, we feel that more detailed simulations of telescope/instrument 
performance and exoEarth signal extraction are required to justify any further reductions. 
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We currently assume one year of total exposure time out of an assumed five year mission. 
A coronagraph could require one additional year of overheads to achieve the desired contrast 
on every target. Thus, our assumptions could implicitly require two years of a five year 
mission. Given that yield weakly depends on total exposure time (0 = 0.33), increases to 
the total exposure time provide little additional yield and would require a larger fraction of 
the total mission. We therefore find it unlikely that increasing the time devoted to exoEarth 
science for a general astrophysics missions could increase yield by more than ~ 10%. 

The weak total exposure time dependence also argues against a multi-telescope coron¬ 
agraph mission, because flying N telescopes is equivalent to one telescope lasting N times 
as long. For example, to double the exoEarth candidate yield, one can either increase the 
aperture by 45% or fly eight identical telescopes. Cost scalings likely make the former a 
more attractive option. 

Under the baseline assumptions of zero detector noise, spectral characterization does 
not strongly impact the yield. This assumes that Rc = 50 spectra are obtained only for the 
true exoEarth candidates and equivalent spectra are not obtained for potentially confusing 
crescent phase gas giants or other background objects. Depending on the observing strategy 
implemented, this may be a reasonable approximation. However, the weak impact of spectral 
characterization time is also a result of the weak dependence of yield on total exposure time. 
Even if half of the total exposure time is devoted to spectral characterization, the mission 
yield is only reduced by 20%. 

If detector noise is non-zero, the increase in the required exposure times can significantly 
impact yield. To define acceptable detector noise parameters, we must choose an acceptable 
level of yield reduction from Figure [101 We choose a 10% reduction in yield as a maximally 
acceptable value, which could be offset by a 5% increase in aperture size. The impact of 
noise from the detection and characterization coronagraphs’ detectors will compound, so we 
allow for at most a 5% reduction in yield due to each of the coronagraphs’ detectors. These 
“allowable” yield reductions are shown as solid lines in Figure [TOl Any noise properties to 
the lower-left of these curves would be acceptable from a yield perspective for our baseline 
mission parameters. 

Because an IFS spreads the planet’s light over a much larger number of pixels, tolerable 
noise properties are much more stringent. For planet detection with an IFS under the 
baseline assumptions, we require a detector with dark current < 0.0001 pix“^ s“^ and read 
noise < 0.2 pix“^. Energy-resolving detectors, which spread the planet’s light over a small 
number of pixels, must only achieve dark currents < 0.01 pix“^ s“^ and read noise < 2 pix“^, 
on par with current CCD noise performance. 
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We note, however, that we have not considered exposnre time limitations dne to pointing 
reqnirements or the motion of the planet being observed. The allowable noise parameters for 
the characterization coronagraph’s detector may become more constrained as one limits the 
spectral characterization time. Further, the allowable noise properties of the characterization 
coronagraph’s detector are sensitive to the values of Rc and S/Nc, as shown by Figure [TT] 
Thus the tolerable detector noise properties should be defined in terms of the most difficult 
spectral line one wants to detect. For our baseline mission in this work, we considered only 
the detection of water near 1 pm, which requires Rc ~ 50 and S/N^ ~ 5. If one assumes the 
mission must be capable of detecting O 2 near 0.75 pm, such that Rc ~ 150 and S/Nc ~ 6 
flBrandt fc Spiegel! l2014j) . then the dot-dashed curves in Figure fTTl should be considered in 
defining the tolerable noise parameters. 


5.1. Lower limits on aperture size 


For our baseline mission, we assume noiseless detectors, relatively high total throughput 
of 20% for a coronagraphic mission, raw detection contrast of 10“^°, and a 360° detection 
zone between the IWA and OWA. Additional sources of performance degradation, like jitter, 
are not considered and implicitly assumed to be negligible. We can therefore consider our 
baseline coronagraph performance optimistic. 

Let’s assume for the moment that such an optimistic mission is attainable and examine 
the yield as a function of telescope aperture under different astrophysical assumptions. In a 
pessimstic/low risk scenario, we would design a mission robust to a low value of p®. Thus, 
for our pessimistic case we choose = 0.1, our Icr lower limit from Table [2J It is extremely 
unlikely that a mission would be planned using anything greater than the expected value for 
? 70 , so for the optimistic/high risk scenario we adopt = 0.16, our expected value for the 
OKHZ. 


Figure dn] shows the exoEarth candidate yield as a function of aperture size and median 
exozodi level under these two assumptions. The current constraint on the median exozodi 


level is < 60 zodis (iMennesson et al.l 1201411 . the rightmost curve in each plot. If LBTI can 
further constrain the median exozodi level, perhaps as low as ~3 zodis, the yield would 
follow the curves on the left. 
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nexoz„d.s = 3 5 10 20 60 n_ = 3 5 10 20 



Fig. 16.— ExoEarth candidate yield vs aperture for the optimistic (higher risk) and pes¬ 
simistic (lower risk) values of i]®. Our minimum acceptable yield for an exoEarth charac¬ 
terizing mission is shown as a dotted line. To exceed the minimum yield under the most 
optimistic (and potentially unrealistic) assumptions, we require an aperture > 8.5 m. 


To place a lower limit on t he required apertu re size, we must choose a minimum ac¬ 
ceptable yield. Equation 3 from IStark et al.l (120141 ) provides a quantitative means to dehne 
a yield goal. Here, we present a modihed form that accounts for the counting statistics of 
all exoEarth candidates, not just the “true” exoEarths: 


^Ec — he 


log(l-g) 

log(l - he/x)’ 


( 8 ) 


where IVec is the target yield of exoEarth candidates, C is the conhdence of Ending at least 
one candidate with x detectable in its atmosphere, fx is the fraction of exoEarth candidates 
that have x detectable in their atmospheres, and x can be water, oxygen, methane, or 
anything else scientifically compelling. 


Stark et ahl (120141 1 suggested designing a mission robust to = 0.1 with 99.7% confi¬ 


dence. To place a lower limit on the aperture size, we relax this constraint to fx = 0.1 with 
95% confidence, which we consider the minimum acceptable yield. Assuming that detectable 
levels of X define an “Earth-like” planet, designing a mission to this standard guarantees, 
with 95% confidence, the detection of at least one Earth-like planet if one in ten Earth-sized 
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planets in the HZ are Earth-like. Snbstitnting these valnes into Eqnation[8]gives a yield goal 
of 30 exoEarth candidates for both the optimistic and pessimistic valnes of r]/^. 

We plot the minimnm scientifically compelling yield goal as a horizontal dashed line in 
Fignre [161 Even in the most optimistic scenario possible, with an optimistic coronagraph, 
optimistic astrophysical assnmptions, and an optimistic observation plan in which spectral 
characterization is only performed on those planets that are trnly exoEarth candidates, we 
reqnire an apertnre > 8.5 m to achieve the minimnm yield goal. Given that the planned 
Space Lannch System (SLS) Block-2 with 10 m fairing, the most capable space lannch 
rocket for the foreseeable fntnre, conld not lannch a monolithic apertnre > 8 m in diameter, 
a coronagraphic exoEarth detecting mission may reqnire a segmented apertnre. 


We note that Fignre [in] assnmes no change in telescope or coronagraph performance 
as a fnnction of D. In reality, segmented apertnres reqnire specialized coronagraphs, so we 
expect some change in performance at the transition from monolithic to segmented apertnres, 
which will occnr at < 8 m given the largest SLS fairing diameter. To date, the majority of 
coronagraph design work has focnsed on monolithic, off-axis apertnres. Cnrrent coronagraphs 
designed to work with segmented apertnres have comparatively lower perfor mance, thongh 


there are designs that, in principle, can perform on-par with off-axis designs flN’Diave et ah 


20151: Ihvonl |2015|) . Thns, if one were to plot yield as a fnnction of apertnre for today’s 


coronagraphs, a significant rednction in yield wonld occnr when crossing the 8 m threshold, 
which wonld eventnally be offset at larger apertnres. This degradation in yield at the 8 m 
transition is not well qnantihed and depends on the parameters of the specihc coronagraphs 
chosen for monolithic and segmented apertnres. We plan to address this issne and the yield 
of specific coronagraphs in fntnre work. 


Finally, we note that improving the constraint on the median exozodi level can snbstan- 
tially rednce the reqnired apertnre. A mission robnst to 60 zodis and = 0.1 reqnires an 
apertnre > 16 m. If LBTI Ends that the median exozodi level < 3 zodis, we wonld reqnire 
an apertnre of ~ 11 m. 


6. Conclusions 

We have developed a mnlti-wavelength completeness calcnlator that handles mnltiple 
visits to every star and simnltaneonsly optimizes the nnmber of visits, delay time between 
visits, and exposnre time of every observation to maximize exoEarth candidate yield for a 
coronagraphic direct-imaging mission. We have npdated onr astrophysical assnmptions based 
on recent literatnre, nsing circnlar orbits within the optimistic habitable zone as defined 
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by iKopparapu et al.l (120131) . For a population of rocky planets residing in this classical 
habitable zone with median radius = 1 i?®, we find p® = 0.16 ± 0.06 by scaling the results 


of lPetigura et al.l fl2013al) . 


We find that optimized revisits increase yields by 35-75%, providin g a ~ 15% reductio n 


in the required aperture size. We find scaling results similar to those of IStark et ahl (1201411 . 
Yield is most strongly dependent on aperture size, scaling roughly as oc D^. We find a 
moderately weak dependence of yield on total exposure time (oc which argues against 

significant increases in exoEarth science time as well as mission concepts with multiple small 
telescopes. 


As long as post-processing enables the reliable detection of point sources brighter than 
26.5 mags dimmer than the host star, the yield is a weak function of the coronagraph’s raw 
contrast and the contrast floor; a factor of 10 degradation in the detection coronagraph’s 
contrast from the baseline value of 10“^*^ reduces yield by ~ 40%. The yield is strongly 
dependent on the inner working angle of the coronagraph (oc IWA“°'®®); a degradation of the 
inner working angle from the baseline value of 2\c/D to 3Xc/D reduces yield by ~ 35%. We 
find a weak to moderately weak dependence on all other high level coronagraph parameters: 
outer working angle, throughput, bandwidth, and PSF sharpness. As long as spectra are 
obtained for a small subset of the observations (e.g., only for likely exoEarth candidates), 
the yield is a very weak function of the characterization coronagraph’s contrast. 


Yield is a weak function of the median exozodi level, but there is considerable room for 
improvement in median exozodi level constraints. Constraining the median exozodi level to 
< 3 zodis, which may be possible with LBTI, could more than double the expected yield 
and reduce the required aperture size by 40%. We find little benefit in knowing the exozodi 
levels of individual stars a priori and find that the precise distribution of exozodi negligibly 
impacts yield for apertures > 4 m. 


We have included the effects of detector noise and spectral characterization time and 
calculated the yield under the assumption that imaging and spectral characterization are 
drawn from the same total time budget. In the limit of noiseless detectors, the required 
spectral characterization time does not significantly impact mission yield. 


If spectral information is not required for wavefront control and a simple broadband 
imager can be used for detection, or if an energy-resolving detector can be used, constraints 
on allowable detector noise are modest, with dark current < 0.01 s“^ pixel“^ and read noise 
;< 2 pixel“^. If an energy-resolving detector is not an option and an IFS is required for 
detection, we require improved detector noise performance, with dark current < 0.0001 s“^ 
pixel“^ and read noise < 0.2 pixel”^. 
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We place lower limits on the aperture size for a number of astrophysical assumptions, 
while assuming optimistic telescope and coronagraph parameters. We conclude that a coro- 
nagraphic direct-imaging mission requires a > 8.5 m aperture to provide a reasonable chance 
of hnding an Earth-like planet. 

Future yield estimates could beneht from an improved target list, vetted in more detail 
for completeness, binarity, and stellar diameter. A more realistic simulation of coronagraph 
performance, including contrast as a function of separation and Airy throughput as a function 
of separation, would enable us to compare different coronagraph designs on an equal footing. 
Finally, future efforts should also focus on analogous yield maximization methods optimized 
for starshade-based mission concepts. 
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